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CHAPTER  I 


INTRODUCTION 

A manifestation  of  the  rocket  exhaust  of  a missile  in  flight  is  the  pre- 
sence of  a trailing  plume.  And  whereas  the  missile  skin  is  made  of  conduct- 
ing metal,  the  rocket  exhaust  plume  consists  of  partially  ionized  gases  or 
plasma  which  electrically  gives  the  plume  region  a finite  conductivity.  When 
this  composite  conducting/partially-conducting  object  is  immersed  in  an  elec- 
tromagnetic environment,  it  is  necessary  to  estimate  the  skin  current  on  the 
missile  with  plume  in  order  to  determine  how  on  board  subsystems  may  respond 
Ll,  2,  3,  4].  In  a laboratory  situation  (e.g.,  in  an  anechoic  chamber), 
however,  it  is  difficult  to  accommodate  the  presence  of  the  plume  because 
of  the  attendant  hazards  and  expenses.  Thus  testing  is  often  done  without 
accounting  for  the  plume.  The  problem  we  address  in  this  report  is  to  esti- 
mate what  differences  one  might  expect  in  subsystem  response,  under  electro- 
magnetic compatibility  (EMC)  testing,  with  and  without  a trailing  plume 
present.  The  consideration  of  the  actual  internal  coupling  mechanism  is  an 
extremely  difficult  problem  whose  consideration  would  take  us  far  afield. 
However,  if  the  coupling  is  primarily  through  electrically  small  apertures  in 
the  missile  skin,  then  it  is  known  that  the  coupling  to  the  interior  is  pro- 
portional to  a linear  combination  of  the  current  and  charge  that  exist  at 
the  aperture  when  it  is  covered  over  by  a metallic  conductor.  Thus,  computa- 
tion of  this  so-called  "short  circuit  current"  on  the  missle  with  and  without 
the  plume  present  should  be  sufficient  to  determine  what  differences  in 
the  coupling  might  be  attributable  to  the  presence  of  the  plume.  The  "short- 
circuit  charge,"  of  course,  can  be  found  from  the  current  using  the  equation 
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of  continuity. 


In  the  past,  thin  wire  theory  has  been  used  to  estimate  the  current 
induced  on  a missile  with  plume  [1-4].  There  the  missile  is  viewed  as  a thin 
conducting  wire  of  radius  'a'  and  length  'h',  as  shown  in  Fig.  1,  and  the 
plume  is  modeled  as  an  attached  lossy  wire  of  the  same  radius  as  the  missile. 
The  radius  a is  assumed  to  be  small  compared  to  the  wavelength  A and  to  the 
axial  length  of  the  missile  and  plume,  thus  permitting  the  use  of  thin  wire 
theory.  The  conductivity  CT  of  the  plume  may  be  constant  or  axially  varying 
depending  on  whether  a homogeneous  or  axially  inhomogeneous  model  is  assumed, 
respectively.  Recent  calculations  using  the  Aero  Chem  low  altitude  plume 
program  (LAPP)  [5]  indicate,  however,  that  the  plume  conductivity  a is 
rather  strongly  inhomogeneous  both  axially  and  radially,  and  that  the  plume 
radius  is  considerably  different  from  that  of  the  missile  and  varies  along 
the  plume  axis.  In  addition,  these  effects  are  also  strongly  influenced  by 
the  missile  altitude.  Therefore,  more  sophisticated  models  and  methods  seem 
to  be  necessary  to  study  the  current  induced  on  the  missile  with  and  without 
the  p^ume. 

In  this  report,  the  thin  wire  model  and  body  of  revolution  models  are 
used  to  assess  the  currents  induced  on  a missile  with  and  without  a trailing 
plume  present.  Chapter  II  summarizes  the  simple  and  efficient  thin  wire 
theory  which  is  useful  for  estimating  order  of  magnitude  effects.  Although 
the  plume  radius  is  fixed  in  this  model  it  is  capable  of  modeling  either  a 
homogeneous  plume  or  an  axially  inhomogeneous  plume.  In  Chapter  III,  several 
body  of  revolution  models  are  described,  compared,  and  used  to  analyze  large 
radius  and  inhomogeneous  plume  models.  Extensive  results  are  ’iven  in  both 


f 


2 


l 


I 


chapters  to  permit  one  to  compare  the  differences  in  currents  induced  on  the 
missile  with  and  without  the  plume  as  functions  of  frequency,  incident  angle 
of  the  plane  wave  source,  plume  conductivity,  and  plume  shape.  In  the  body 
of  revolution  models,  only  the  circumferentially  uniform  Fourier  component 
of  the  current  has  been  computed.  For  reasonably  thin  missiles  the  non- 
symmetric  current  modes  are  only  weakly  coupled  to  the  plume.  Hence,  these 
modes  of  missile  current,  although  sometimes  significant  in  themselves,  are 
not  expected  to  change  much  with  or  without  the  plume  present. 


CHAPTER  II 
THIN  WIRE  MODEL 

In  this  chapter,  a numerical  determination  of  the  induced  currents  on  a 
missile  with  a trailing  plume  is  made  using  simple  wire  models.  In  the  analy- 
sis, the  missile  with  plume  is  viewed  as  a thin  conducting  wire  of  radius 
a <<  A,  as  shown  in  Fig.  1,  with  a resistively  loaded  section,  thus  permitting 
the  use  of  thin  wire  theory  [6,7,8].  Both  homogeneous  and  axially  inhomoge- 
neous plumes  are  considered.  Numerical  results  for  the  induced  currents  are 
given  as  a function  of  frequency,  incidence  angle  of  the  plane  wave  source, 
and  the  plume  conductivity.  The  main  objective  is  to  examine  the  differences 
in  currents  induced  on  a missile  with  or  without  the  plume.  The  presence 
of  canards  and  fins  might  be  accounted  for  by  introducing  attached  cross- 
wires into  the  model.  However,  a typical  result  for  a 10'  missile  with  8" 
crossed  wires  is  as  shown  in  Fig.  2,  where,  as  can  be  seen,  other  than  the 
expected  discontinuity  in  current  at  the  junction,  the  current  distribution 
and  peak  currents  are  not  significantly  changed.  Thus  the  effects  of  body 
appendages  are  considered  to  be  minor,  and  will  not  be  included  throughout 
the  remainder  of  these  studies. 

2.1.  Theory 

Consider  the  thin  wire  model  as  depicted  in  Fig.  1.  The  missile  (with 
radius  "a"  and  length  ”h")  is  regarded  as  a conducting  wire  and  the  plume 
(with  radius  "a",  length  and  conductivity  "a")  is  regarded  as  an  impe- 

dance (or  resistively)  loaded  section.  Assuming  the  current  flows  only  in 
the  axial  or  z direction,  one  may  obtain  the  following  integro-dif ferential 
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missile  I plume 


Fig.  1.  Thin  wire  modeling  of  a missile  with  plume. 
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equation  [6,7 ] via  enforcing  the  boundary  conditions  on  the  wire: 


k2  \ dz2  °l 


A (z)  + E1(z)  = Et0tal  = Z(z)I(z) 
Z Z z 


(1) 


where  Z(z)  is  the  equivalent  impedance  along  the  wire  and  is  given  by 

0 , ....  for  the  conducting  region, 

Z(z)  = <!  , , _ . J (Ta) 


(2) 


— T" ( "T  ) ~T°/~  \ , •••  for  the  lossy  region, 

ira  a(z)  ' L ' J1U  ' 


and 

x = d-j)  j - o) 

where  6 is  the  skin  depth  of  the  material.  The  magnetic  vector  potential 
A^Cz)  is  given  by 


Az(z)  =47  J I (z ' )K(z-z ' )dz ' , (4) 

wire 

K(C)  = exp  (-jkQR)/R  (5) 

R = U2  + aZ)H  (6) 

I(z)  ip  the  unknown  axial  current  to  be  determined,  and  E1(z)  is  the  known  in- 

z 

cident  electric  field.  To  solve  Eq.  (1)  the  method  of  moments  is  employed. 
First,  Eq.  (1)  is  tested  with  a set  of  triangle  functions,  then  the  unknown 
current  I(z)  is  expanded  in  a set  of  flat  pulses  with  unknown  coefficients 
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I n = 1,  2,  ...  N.  This  results  in  a set  of  N independent  algebraic  equa- 
n 

tions.  The  unknown  current  coefficients  1^  are  then  found  using  standard 
matrix  inversion  methods.  The  details  are  given  in  references  [7,8]  and 
will  not  be  repeated  here. 

2.2.  Numerical  Results 

Based  on  the  analysis  in  the  previous  section,  an  efficient  FORTRAN  IV 
computer  program  was  developed  to  find  the  unknown  axial  current  on  the  mis- 
sile and  plume..  Fig.  3 shows  the  wire  code  results  for  the  real  and  imagi- 
nary part  of  the  currents  on  a missile  (radius  a.=  .02A,  length  h = 1.11A) 
with  plume  (of  same  radius  and  length  as  the  missile,  e = 1,  o ■ 0.1  u/m) 
illuminated  by  a broadside  incident  plane  wave  (freq.  = 19.6  MHz).  Here  the 
wire  code  results  agree  very  well  with  the  numerical  results  computed  by 
Taylor  [1].  Once  the  validity  of  this  wire  code  has  been  established,  one 
may  readily  study  the  effects  on  the  missile  currents  due  to  the  presence  of 
the  plume  when  the  frequency,  plume  conductivity,  and  incident  angle  of  the 
plane  wave  sources  are  varied. 

A three-dimensional  view  of  the  current  amplitude  on  a missile  (h  = 10', 
a = 2.5")  and  plume  (£  = 20')  is  plotted  in  Fig.  6 as  functions  of  the  plume 
conductivity  (a  = 0,  0.01,  0.05,  0.2,  0.5,  1 and  2 v/m)  and  the  distance 
along  the  axis  of  the  missile-plume  configuration  with  a broadside  incident 
plane  wave  (freq.  = 50  MHz).  The  case  of  zero  plume  conductivity  represents 
a missile  with  no  plume  present.  One  may  notice  that  when  0 is  0.01  zr/m 
there  is  little  effect  caused  by  the  presence  of  the  plume.  When  0 is  in- 
creased to  1 or  2 v/m,  however,  then  rather  significant  differences  in  cur- 
rents take  place.  Plots  for  other  frequencies  (e.g.,  20,  35,  60,  70,  90,  100 
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and  120  MHz)  are  shown  in  Figs.  4-5  and  7-12,  for  a broadside  incidence  plane 
wave.  For  an  incidence  angle  (8^)  of  120°,  Figs.  13-15  show  the  current 
amplitudes  at  20,  50,  and  100  MHz,  respectively,  as  functions  of  the  plume 
conductivity  and  positions  along  the  axis  of  the  missile  and  plume.  Figs. 
16-19  present  current  data  for  the  case  of  150°  incidence  angle  at  35,  50, 
100,  and  120  MHz.  Figs.  20-42  show  current  plots  for  other  angles  of  inci- 
dence at  20,  35,  50,  100,  and  120  MHz. 

The  frequencies  considered  may  be  divided  into  two  categories;  those 
below  the  first  missile  resonance  and  those  above  the  first  resonance,  where 
the  first  missile  resonance  occurs  at  about  50  MHz.  For  frequencies  below 
resonance,  increasing  the  plume  conductivity  essentially  increases  the  elec- 
trical length  of  the  missile  and  hence  invariably  results  in  larger  missile 
currents.  At  resonance  and  above,  however,  increasing  the  plume  conductivity 
typically  "detunes"  the  missile  resonance  currents  and  reduces  the  current. 
With  further  increases  in  conductivity,  however,  the  missile-plume  configura- 
tion may  be  nearly  resonant  in  a higher  order  resonance,  in  which  case  the 
missile-plume  currents  again  increase.  For  incidence  angles  0 of  about 
150°  or  greater,  however,  any  decreases  in  current  are  negligible  so  that  the 
effect  of  increasing  conductivity  is  essentially  to  increase  the  missile 
current.  Many  of  these  phenomena  have  been  reported  in  other  independent 
studies  [4]. 

The  effect  of  the  incident  angle  0^  on  the  missile  current  distribution 
may  also  be  plotted  in  a polar  diagram.  Figs.  43-45  show  the  current  magni- 
tude at  the  middle  of  the  missile  with  and  without  plume  for  frequencies  of 
50,  100,  and  120  MHz,  respectively.  The  current  magnitudes  at  3'  from  the 


Fig.  4.  Contour  plot  of  current  magnitude  as  functions  of  the  length  along 


the  missile-plume  axis  and  the  plume  conductivity  a (Broadside 
incidence,  f = 20  MHz). 
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Fig.  45.  1 1 1 at  middle  on  the  missile  (f  = 120  MHz). 
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nose  of  the  missile  are  shown  in  Figs.  46-48.  As  can  be  seen,  significant 
changes  in  current  magnitude  occur  as  the  plume  conductivity  increases. 

These  plots  are  also  useful  for  assessing  the  coupling  effect  through  aper- 
tures in  the  missile  skin. 

In  summary,  the  effect  of  a homogeneous  plume  on  the  missile  current 
distribution  has  been  studied  in  this  section  using  the  thin-wire  modeling 
and  analysis.  Extensive  data  show  that  minimal  effects  are  observed  if  the 
plume  conductivity  a is  less  than  0.01  v/m,  but  if  a > 1 u/m,  then  signifi- 
cant differences  are  present  as  compared  to  the  case  of  no  plume  present. 

If  the  missile  is  short  with  respect  to  the  wavelength  (e.g.,  for  a 10'  mis- 
sile at  20  MHz)  the  plume  usually  enhances  the  current  magnitude  on  the  mis- 
sile. However,  if  the  missile  is  near  or  above  resonance  (e.g.,  for  a 10' 
missile  at  frequencies  larger  than  35  MHz) , the  plume  usually  suppresses  mis- 
sile current  for  incident  angles  0^  <.  120°,  while  for  8^  >.  150°  enhancement 
of  current  occurs  again. 

In  addition  to  the  above  results,  the  thin  wire  code  may  also  be  used  to 

assess  the  effect  of  an  axially  inhomogeneous  plume  with  an  axially  varying 

conductivity  as  shown  in  Fig.  49.  Note,  however,  that  in  the  calculation 

of  the  equivalent  impedance,  the  radius  of  the  plume  may  differ  from  that  of 

the  missile,  but  in  the  analysis  we  assume  the  plume  to  have  the  same  radius 

as  the  missile.  Figs.  50  and  51  show  the  current  magnitude  on  the  missile 

with  and  without  the  inhomogeneous  plume  at  f = 100  MHz  for  0^  = 90°  and 

165°  respectively.  As  in  the  homogeneous  plume  case,  minimal  effects  are 

present  if  the  maximum  plume  conductivity  a <.0.01  v/m.  For  higher  plume 

iqbx 

conductivities  and  for  wider  plumes,  more  sophisticated  models  and  analyses 
are  needed  as  discussed  in  the  following  chapter. 
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Fig.  50.  Current  amplitude  |l|  on  a missile  (h  =*  10',  a = 2.5')  with  and 
without  the  inhomogeneous  plume  of  Fig.  49  (broadside  incidence, 
f = 100  MHz). 


0 10  20  30 

z (ft) 


Fig.  51.  | I | on  a missile  with  and  without  the  inhomogeneous  plume  of 

Fig.  49  (6i  = 165°,  f - 100  MHz). 
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CHAPTER  III 


BODY  OF  REVOLUTION  MODEL 


The  plume  from  a rocket  exhaust  tends  to  broaden  away  from  the  missile 
axis,  particularly  at  higher  altitudes,  so  that  thin-wire  plume  models  are 
often  inadequate.  Hence  we  consider  in  this  chapter  a more  sophisticated 
missile/plume  model.  This  model  consists  of  a conducting  cylinder  modeling 
the  missile,  with  a homogeneous  or  inhomogeneous  lossy  dielectric  body  of 
revolution  trailing  behind.  For  this  composite  body  of  revolution  scattering 
problem  [9],  the  skin  currents  on  the  missile  and  plume  may  be  calculated 
using  either  a Surface  Integral  Equation  (SIE)  [ 10]  or  a so-called  Approxi- 
mate Integral  Equation  (AIE)  [ 18  ] technique.  The  analysis  for  a homogeneous 
plume  model  involves  a moment  method  solution  of  a set  of  coupled  vector 
integral  equations  and  is  summarized  in  Section  3.1.  Since  the  radius  of 
the  missile  is  typically  small  in  terms  of  wavelengths,  the  currents  on  the 
missile  consist  primarily  of  the  n=0  and  n«l  circumferential  Fourier  modes. 
Although  either  or  both  of  these  modes  can  make  significant  contributions  to 
the  current  at  various  points  on  the  missile,  it  is  probable  that  only  the 
n=0  mode,  i.e.,  the  circumferentially  uniform  mode,  is  tightly  coupled  to  the 
plume.  Since  we  are  interested  primarily  in  the  disturbance  of  the  missile 
current  due  to  the  presence  of  the  plume,  and  this  is  likely  to  have  no  ef- 
fect on  the  n=l  mode,  we  consider  here  only  the  circumferentially  uniform  mode. 
In  Section  3.2  the  scattering  of  an  inhomogeneous  plume  model  is  considered 
using  the  Approximate  Integral  Equation  (AIE)  approach.  Numerical  results 
obtained  using  these  methods  are  presented  in  Section  3.3.  The  accuracy  and 
applicability  of  the  AIE  method  is  also  discussed. 
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. Surface  Integral  Equation  (SIE)  Method 


Consider  the  composite  metal-lossy  dielectric  body  of  revolution  (BOR) 
model  depicted  in  Fig.  52.  The  conducting  portion  represents  the  missile 
and  the  dielectric  portion  represents  the  homogeneous  plume.  From  Maxwell's 
Equations,  Green's  theorem,  and  the  boundary  conditions  [10,11]  one  readily 
obtains  the  following  coupled  vector  integral  equations  for  the  base-contact 
BOR(I)  model: 
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1 J o 1 
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The  above  integral  equations  are  implicitly  valid  for  tangential  components 
only,  and  may  also  be  derived  via  the  surface  equivalence  principle.  3 

1 

and  3_  are  the  unknown  electric  current  densities  on  the  air-conductor 
S2 

interface  S,  and  the  conductor-dielectric  interface  S„,  respectively.  J 

1 ^ o 
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and  K are  the  unknown  electric  and  magnetic  current  densities  on  the  air- 

S3 

dielectric  interface  S . The  boundary  conditions  also  require  that  K 

j 

vanishes  and  J , J and  J are  equal  at  the  junction  of  the  air,  conductor 

1 ,2  _3  _ 

and  dielectric.  E (r)  and  H1(r)are  the  incident  electric  and  magnetic  fields 
and  <J>2  = exp(-jk2  | r-r ' | )/ j r-r ' | and  = exp(-jkQ |r-r' |)/|r-r' | are  the 
Green's  functions  in  the  lossy  dielectric  and  in  free  space,  respectively. 

The  wavenumbers  are  k„  = oj/y  e.  , k = oo/y  £ , and  e-  is  the  (complex) 

L Oil  o oo  i. 

permittivity  of  the  dielectric.  The  time  dependence  is  assumed  to  be  exp(ju)t). 


For  the  point-contact  BOR(II)  model,  one  may  readily  formulate  the 
following  integral  equations: 
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(13) 


One  notes  that  for  this  model  the  conductor-dielectric  interface  reduces 

to  a point.  Thus  J vanishes  and  J and  J both  approach  zero  at  the 

2 1 3 

contact  point.  All  other  quantities  are  defined  the  same  as  in  Eqs.  (7) 

(10  ). 

In  order  to  evaluate  the  unknown  currents, one  first  expands  the  currents 
in  a Fourier  series  in  the  circumferential  angle  <j>  and  uses  triangle-functions 
to  represent  the  variation  of  the  current  with  the  coordinate  t.  The  unknown 
coefficients  of  the  triangular  expansion  functions  are  then  determined  from 
a set  of  algebraic  equations  which  are  obtained  by  testing  the  above  equa- 
tions ( 7 ) - (13)  with  appropriate  triangle  testing  functions  and  equating 
the  (^'independent  Fourier  components  on  both  sides  of  the  equations.  The 
solution  procedure  is  rather  involved  but  essentially  parallels  that  found 
in  [10,  12],  except  that  half-triangles  are  used  to  represent  the  nonzero 

currents  J , J , and  J at  the  air- conductor-dielectric  interface  in  the 
S1  2 3 

BOR(I)  model.  The  above  SIE  method  has  been  successfully  applied  to  the 
composite  body  represented  by  the  missile  with  a homogeneous  plume  model 
and  the  numerical  results  are  shown  in  Section  3.3.  In  principle,  the  SIE 
method  can  also  be  applied  to  the  composite  body  consisting  of  a missile 
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with  a multi-layered  inhomogeneous  plume  by  using  an  invariant  imbedding 
technique  [13,14].  This  approach  is  currently  under  investigation. 

3.2.  Approximate  Integral  Equation  (AIE)  Method 

It  is  a very  difficult  and  complicated  matter  to  accurately  model  a 

missile  with  an  attached  plume,  whose  composition  varies  both  axially  and 

radially.  However,  an  approximate  so-called  AIE  method  is  employed  here  to 

obtain  a simple  first  order  solution  to  this  problem.  The  approach  here  is 

to  derive  an  approximate  boundary  condition  relating  the. tangential  fields 

at  the  outer  surface  of  the  plume  and  to  model  the  effect  of  the  plume  by 

an  impedance  type  boundary  condition  on  the  fields  at  the  plume  surface. 

Specifically,  the  inhomogeneous  plume  is  first  divided  into  N segments  along 

its  axis.  To  estimate  the  equivalent  surface  impedance  at  each  segment,  one 

then  calculates  the  appropriate  ratio  of  surface  electric  and  magnetic  fields. 

I * i 

To  obtain  these  fields,  we  locally  approximate  the  radially  inhomogeneous 

segment  by  a multilayered  infinite  circular  cylinder,  which  is  then  solved 
via  an  iteration  procedure  developed  by  Bussey  and  Richmond  [15].  Once  the 
surface  impedance  is  determined,  the  unknown  electric  and  magnetic  current 
densities  on  the  outermost  surface  of  the  plume  are  thus  related  to  each 
other  through  this  surface  impedance  or  so-called  Leontovich  boundary  condi- 
tion [16,17], 

Ks  x n = ZSJS  (14) 

where  Z is  the  calculated  surface  impedance  of  the  inhomogeneous  plume  seg- 

D 

ment,  and  n is  the  unit  normal  vector  to  the  plume  surface. 

The  use  of  the  Leontovich  or  surface  impedance  approximation  relies  on 
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the  following  assumptions:  (a)  the  ratio  of  surface  electric  and  magnetic 
fields  is  not  a strong  function  of  the  angle  of  incidence,  (b)  there  is  not 
a strong  interaction  between  adjacent  subsections  and  (c)  the  plume  surface 
may  be  approximated  locally  as  a cylinder.  Conditions  (a)  and  (b)  should 
be  satisfied  for  moderate  and  high  conductivities  near  the  plume  surface. 
Condition  (c)  should  hold  everywhere  except  possibly  in  the  missile-plume 
junction  region. 

Since  the  impedance  condition  ( 14  ) relates  the  surface  fields  or  equi- 
valent currents,  the  number  of  unknown  quantities  is  reduced  and  one  is  able 
to  derive  an  AIE  for  only  one  of  the  unknown  plume  surface  currents.  Similar 
to  the  derivation  of  Eqs.  ( 7 - 10  ) , with  condition  ( 14  ) , one  may  easily 
arrive  at  the  following  AIE: 
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(15) 


The  above  equation  is  implicitly  valid  for  tangential  components  only.  Note 
that  Zg  is  zero  on  the  conducting  rocket  and  varies  with  position  along  the 
outer  surface  of  the  plume.  The  integral  equation  may  be  solved  using  the 
same  numerical  procedures  indicated  in  the  previous  section.  Here  again,  only 
the  zeroth  Fourier  component  of  the  unknown  current  Jg  is  needed  on  the  mis- 
sile structure. 

The  AIE  method  may  also  be  applied  to  homogeneous  plume  models  where  the 
surface  impedance  is  calculated  according  to  Eq.  ( 2 ).  Numerical  results 
for  homogeneous  plume  models  may  thus  be  obtained  using  either  the  SIE  or 
the  AIE  method,  and  some  results  using  both  methods  are  given  in  the  next 
section  for  comparison.  The  applicability  and  accuracy  of  AIE  method  are 
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also  discussed. 

3.3.  Numerical  Results 

The  SIE  solution  for  the  total  axial  current  on  a homogeneous  dielectric 
circular  cylinder  (length  h = . 737A,  radius  a » .025A,  £r  = 50,  conductivity 
a = 20.37  y/m)  is  plotted  in  Fig.  5.3  for  the  case  of  a broadside  incident 
plane  wave  (f  = 9.45GHz).  The  good  agreement  among  the  various  numerical 
solutions  validates  the  calculation  of  the  zeroth  order  Fourier  component 
using  the  SIE  method.  In  Figs.  54  and  55  the  current'  magnitudes  on  various 
missile-plume  models  are  plotted  for  broadside  and  30°  incident  plane  wave 
illuminations  (f  = 50  MHz),  respectively.  For  the  base-contact  BOR(I)  model, 
the  SIE  solutions  are  in  good  agreement  with  the  thin-wire  code  solutions. 

The  SIE  solution  for  the  point-contact  BOR(II)  model,  however,  yields  sig- 
nigicantly  different  results  from  the  thin-wire  model  for  the  broadside 
incident  case.  The  difference  is  obviously  due  to  the  fact  that  in  this 
model  current  is  not  allowed  to  flow  from  the  missile  into  the  plume  region 
because  the  contact  area  between  the  missile  and  plume  is  zero.  In  the  30° 
incident  case,  the  current  in  the  junction  region  happens  to  be  small  anyway, 
so  that  requiring  it  to  vanish  at  the  interface  in  the  BOR(II)  model  does  not 
significantly  change  the  overall  current  magnitude.  In  Figs.  56  and  57 
we  show  the  current  magnitude  on  various  missile-plume  models  with  differing 
plume  shapes  for  broadside  and  30°  incident  plane  waves  (f  - 50  MHz),  re- 
spectively. In  each  case  the  plume  is  homogeneous  and  the  conductivity  is 
0.2  TT/m.  As  in  the  case  of  the  thin-wire,  at  broadside  incidence,  the 
plume  suppresses  the  missile  current  amplitudes,  whereas  near  grazing  inci- 
dence, e.g.,  0^  < 30°,  the  current  is  enhanced.  Some  influence  due  to  the 
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Figure  55.  Current  magnitudes  on  the  same  missile  as  Figure  54  but  fo 
angle  of  incidence  0.  = 30°. 


plume  shape  and  length  can  also  be  seen.  In  general,  for  a fixed  conduc- 
tivity, the  larger  the  plume  radius,  the  more  nearly  it  behaves  like  a 
perfect  conductor.  This  is  because  increasing  the  plume  cross  section  inc- 
reases the  region  of  current  flow,  in  effect  reducing  the  impedance  per  unit 
length. 

A comparison  of  the  current  obtained  both  by  the  SIE  and  the  AIE  methods 
for  a missile  with  a homogeneous  plume  (a  = 0.1  and  1 TT /m)  is  made  in  Figs. 

58  and  59  for  a broadside  incident  plane  wave  (f  = 50  MHz  and  100  MHz,  res- 
pectively) . It  is  noted  that  for  higher  plume  conductivity  or  frequency  the 
AIE  solution  should  yield  better  agreement  with  the  SIE  solution,  as  the 
figures  verify.  As  the  plume  conductivity  or  the  frequency  decreases,  how- 
ever, one  may  expect  that  the  simpler  AIE  solution  will  not  be  as  accurate 
as  the  SIE  solution.  Nonetheless,  the  AIE  method  yields  a first  order  solu- 
tion for  the  inhomegeneous  plume  model. 

To  demonstrate  the  applicability  of  AIE  solution  we  first  consider  the 
plume  conductivity  profile  of  a Redeye  rocket  at  5000  ft.,  altitude  as  shown 
in  Fig.  60.  (this  conductivity  profile  has  been  computed  by  means  of  the 
so-called  LAPP  computer  code  at  the  Missile  Research  and  Development  Command, 
Redstone  Arsenal,  Alabama  [5].)  Several  views  of  the  profile  are  presented 
in  order  to  observe  the  detailed  structure.  Fig.  61  shows  the  maximum 
radius  used  in  our  model  of  the  plume  at  each  point.  At  points  on  the  sur- 
face in  Fig.  61  , an  equivalent  surface  impedance  was  computed  and  used  in 

the  AIE.  The  resulting  AIE  current  solutions  for  this  inhomogeneous  plume 
model  are  plotted  in  Figs.  62  and  63  for  the  cases  of  broadside  and  30° 
incident  plane  waves  (f  = 100  MHz),  respectively.  Also  plotted  are  the  AIE 
solutions  for  a homogeneous  plume  model  with  plume  conductivity  O = 0.1  mhos4n. 
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Figure  61.  Maximum  rad 
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Figure  63.  Current  magnitude  on  the  same  configurations  as  Figure  62 
but  with  angle  of  incidence  9^  = 30°. 
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Since  this  conductivity  is  considerably  higher  than  the  average  conductivity 
of  the  inhomogeneous  model,  one  would  expect  the  currents  on  the  plume  in 
the  homogeneous  model  to  be  slightly  higher  than  those  of  the  inhomogeneous 
case,  as  is  confirmed  by  Figs.  62  and  63.  Note  that  the  currents  on  the 
conductor  portion  are  almost  equal,  however.  For  this  example,  then,  the 
homogeneous  model  might  be  adequately  substituted  for  the  actual inhomogenous 
plume.  Indeed,  for  most  of  the  data  we  have  considered  here  the  results  do 
not  seem  to  depend  strongly  on  the  actual  model  employed. 

Consider  next  the  inhomogeneous  plume  corresponding  to  a Chaparral 
missile  (h  = 3m,  a = 6.35cm)  at  a 5000  ft.  altitude,  with  the  con- 
ductivity profile  shown  in  Fig  64  . The  corresponding  maximum  conductivity 
profile  is  plotted  in  Fig.  65  along  the  plume  axis.  One  notes  that  the 
plume  conductivity  of  the  long  plume  model  is  rather  low  near  the  nozzle 
(i.e.  z < 0.684  m) . It  is  agreed  Chat  computed  plume  conductivities  gener- 
ated by  the  LAPP  code  in  the  nozzle  region  are  rather  inaccurate  due  to  an 
inability  to  adequately  model  the  complicated  plume  dynamics  near  the  nozzle. 
In  order  to  assess  the  effect  of  the  inaccurate  modeling  in  the  nozzle  region, 
we  have  adopted  an  alternative  plume  model  in  which  the  low  conductivity 
section  near  the  rocket  nozzle  has  been  removed  (see  Fig.  65  ).  This  model 

is  termed  the  "short  plume  model."  Figs.  66  and  67  show  the  current 
distribution  on  these  two  models  illuminated  by  a plane  wave  at  50  MHz  with 
incident  angle  0^  equal  to  90°  and  30°,  respectively.  As  can  be  seen,  the 
two  inhomogeneous  plume  models  give  rather  different  results.  The  long 
plume  model  with  the  low  conductivity  section  near  the  rocket  nozzle  (at  3m.) 
tends  to  reduce  the  current  amplitude  near  the  missile  tail,  while  the  short 
plume  permits  large  curents  to  flow  on  to  the  plume  section,  in  this  case 
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66.  Current  magnitude  on  a Chaparral  missile  (h  = 3m,  a = 6.35  cm) 
with  and  without  the  two  inhomogeneous  plume  models,  0 = 90°, 

f = 50  MHz. 


reducing  the  peak  current  on  the  missile.  As  has  been  found  for  the  homo- 
geneous plume  models  at  this  frequency,  both  inhomogeneous  plume  models 
show  a larger  current  peak  than  the  no-plume  case  when  the  plane  wave  is  at 
grazing  incidence  (i.e.  for  9^  < 30°).  Figs.  68  to  82  show  plots  of  the 
current  magnitude  along  the  axis  of  the  Chaparral  missile  (h  = 3 m,  £.  = 7.2m, 
a = 6.35  cm)  illminated  by  an  incident  plane  wave  for  a variation  in  fre- 
quency of  30  MHz  to  150  MHz.  The  plots  are  arranged  in  groups  of  three 
showing  results  for  the  (1)  no-plume,  (2)  the  long  plume,  and  (3)  the  short 
plume  case  for  angles  of  incidence  0,  = 30°,  60°,  90°,  120°,  and  150°. 

Again,  the  plots  show  that  because  of  the  low-conductivity  region  at  the 
plume-nozzle  interface  in  the  long  plume  model,  no  really  significant  dif- 
ferences in  missile  current  with  or  without  the  plume  are  observed.  In  the 
short  plume  model,  however,  some  differences  are  noted.  Because  of  the  much 
smaller  plume  conductivities  for  the  Chaparral  rocket,  however,  the  effects 
of  the  plume  are  less  pronounced  than  for  the  Redeye  rocket. 

The  data  indicates  similar  conclusions  for  the  Chaparral  model  as  for 
the  Redeye  model  studied  by  Smith  et.  al  [ h ] , namely,  that  the  missile 
current  is  increased  by  the  presence  of  the  plume  at  lower  frequencies 
because  of  the  increased  effective  length  of  the  nissile^lume  scattering 
configuration.  The  resonance  frequencies  are  slightly  shifted  by  the  plume 
and  some  damping  takes  place  due  to  the  plume  losses.  Note  in  the  figures, 
however,  that  possibly  some  resonance  peaks  have  been  missed  because  of  the 
relatively  large  interval  between  the  frequencies  in  the  data.  Intermediate 
frequency  data  points  were  not  computed  because  of  the  high  cost  of  the 
computations . 
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Figure  74  . Current  on  a Chaparral  missile  as  a function  of  frequency, 
no  plume,  0^  = 90°. 
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Figure  81 


Current  on  a Chaparral  missile  as  a function  of 

long  plume  model,  0 = 150°. 
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CHAPTER  IV 


CONCLUSIONS 


The  effects  of  an  exhaust  plume  on  the  missile  current  distribution  have 
been  studied  via  analyzing  both  thin-wire  and  body-of-revolution  models  of 
the  missile-plume  configuration.  It  was  found  that  the  missile  skin  currents 
are  dependent  on  the  plume  conductivity  and  shape  and  the  frequency  and  the 
incident  angle  of  the  plane  wave  source.  For  the  missile  and  plume  parameters 
considered  here,  minimal  effects  are  observed  if  the  plume  conductivity  'a' 
is  less  than  0.01  mhos/m,  but  if  a > 1 mho/m,  then  significant  differences  in 
the  missile  current  are  produced  as  compared  to  the  case  of  no  plume  present. 
Some  influence  due  to  the  plume  shape  and  length  has  also  been  observed  in 
body-of-revolution  models.  In  general,  for  a fixed  conductivity,  the  larger 
the  plume  radius,  the  more  nearly  it  behaves  like  a perfect  conductor.  This 
is  because  increasing  the  plume  cross  section  increases  the  region  in  which 
the  current  can  flow,  thereby  increasing  the  conductivity  per  unit  length. 

At  resonant  frequencies  and  above,  the  resonant  frequencies  are  slightly 
shifted  by  the  presence  of  the  plume  and  damping  of  the  resonance  peaks  takes 
place  due  to  the  plume  losses.  At  frequencies  below  the  first  missile  reso- 
nance, however,  the  missile  current  is  increased  because  of  the  increased 
effective  length  of  the  missile/plume  scattering  configuration. 

Two  realistic  exhaust  plumes,  with  conductivity  varying  both  axially 
and  radially  in  the  plume  region,  are  used  to  analyze  the  Redeye  and  Chaparral 
missile  plumes  using  the  body-of-revolution  model  and  the  Approximate  Integral 
Equation  (AIE)  method.  The  major  limitation  of  that  analysis  is  that  the 
plume  conductivity  should  be  high  enough  so  the  impedance  condition  (14) 


is  applicable  to  the  plume  surface.  An  accurate  solution  not  relying  on  the 
impedance  condition  approximation  is  currently  under  investigation. 
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